Signatures of spin pairing in a quantum dot in the Coulomb blockade regime 
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Coulomb blockade resonances are measured in a GaAs quantum dot in which both shape defor- 
mations and interactions are small. The parametric evolution of the Coulomb blockade peaks shows 
a pronounced pair correlation in both position and amplitude, which is interpreted as spin pairing. 
As a consequence, the nearest-neighbor distribution of peak spacings can be well approximated 
by a smeared bimodal Wigner surmise, provided that interactions which go beyond the constant 
interaction model are taken into account . 

PACS numbers: 73.20.My, 73.23.Hk, 05.45. +b 
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Recently, the Coulomb blockade (CB) of electronic 
transport through quantum dots, defined in two- 
dimensional electron gases in semiconductor heterostruc- 
tures, has been of considerable interest H. One reason 
is that such dots are model systems to investigate the 
interplay between chaos and electron-electron (e-e) inter- 
actions. Here, a key feature is the distribution of nearest- 
neighbor Coulomb blockade peak spacings (NNS), which 
random matrix theory Q (RMT) predicts to follow a bi- 
modal Wigner surmise P{s) for a non-interacting quan- 
tum dot of chaotic shape, i.e. 
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P(s) = -[S( S )+PP(s)] 



(1) 



P"(s) is the Wigner surmise for the corresponding Gaus- 
sian ensemble, i.e. /3 = 1 for systems with time inversion 
symmetry (Gaussian orthogonal ensemble - GOE), and 
(3 = 2 when time inversion symmetry is broken (Gaussian 
unitary ensemble - GUE). The peak spacing s is measured 
in units of the average spin-degenerate energy level spac- 
ing A = 2ith 2 1 \m* A) , where m* denotes the effective 
mass, and A the dot area. The <5-function in P(s) takes 
the spin degeneracy into account. RMT further predicts 
the standard deviation for P(s) to be a = 0.62 for j3 = 1, 
and a = 0.58 for = 2, respectively [||. 
The comparison to experimental data is made by ap- 
plying the constant-interaction model |^,^| , which allows 
to separate the constant single-electron charging energy 
from the fluctuating energies of the levels inside the dot. 
In disagreement with the predictions of RMT, the exper- 
imentally obtained NNS distributions are usually best 
described by a single Gaussian with enhanced values of 
a |] f|. The data thus look as if spins are absent, al- 
though in Ref. (To), a spin pair has been observed. This 
apparent absence of spins and the different shape of P(s) 
have triggered tremendous recent theoretical work. One 
possible explanation arc additional c-c interactions inside 
the dot ]^, |TT^[T6[ , w hich lead to "scrambling" of the en- 
ergy spectrum |17| , |l0[ ] and can be characterized by the 
interaction parameter r s , defined as the ratio between 



the Coulomb interaction of two electrons at their aver- 
age spacial separation, and the Fermi energy |TT| , ^ -^6[ . 
It is theoretically expected that the NNS distribution be- 
comes Gaussian due to e-e interactions, and that a in- 
creases for r s > 2. 0,^5|. However, all experiments so 
far have been carried out in a regime where an increase in 
a is not expected, i.e. in samples with 0.93 < r s < 1.35 
|B-@|, with the exception of Ref. fl, where r s = 2.1. 
Gate- voltage induced shape deformations of the dot can 
modify the NNS distribution as well. The deformation 
can be described by a parameter x, which corresponds 
to the distance between avoided crossings induced by the 
deformation, measured in units of the CB peak spacing. 
For x ~ 1, the NNS distribution of partly uncorrelated 
energy spectra is measured, resulting again in a Gaussian 
shape with enhanced a fll8 19[. Whether shape defor- 



mations or interactions dominate the shape of the NNS 
distribution is not clear, although there is experimental 
evidence that x < 1 and interactions are more important 

IS 

Here, we report measurements on a quantum dot in which 
shape deformations as well as r s are reduced. We observe 
a pronounced pairwise correlation of both position and 
amplitude of the Coulomb blockade resonances, which is 
sometimes interrupted by kinks in the parametric evolu- 
tion, among other features. We interpret the pairing as 
a spin signature: the energies of two states belonging to 
the same spatial wave function with opposite spin differ 
by an average interaction energy £, which fluctuates with 
a standard deviation of o^, both of which are of the or- 
der of A. We conclude that in previous experiments, spin 
pairing was difficult to observe because it was frequently 
destroyed by avoided level crossings. Furthermore, we 
suggest that the measured NNS distribution can be fit- 
ted to a modified bimodal Wigner surmise, with £ and 
<j£ as fit parameters. 

The sample is a shallow Ga[Al]As heterostructure with 
a two-dimensional electron gas (2DEG) 34 nm below the 
surface. The quantum dot is defined by local oxidation 
with an atomic force microscope p3] (inset in Fig. 1(a)). 
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The lithographic dot area is 280 nm x 280 nm. The dot 
can be tuned by voltages applied to a homogeneous top 
gate and to the planar gates I and II. In order to reduce 
r s as much as possible, we chose a heterostructure with 
a high electron density, further increased by a top gate 
voltage of +100 mV to n e = 5.9 • 10 15 m~ 2 . This results 
in r s = 0.72, which is smaller than in all previous exper- 
iments. Additional screening is provided by the top gate 
pjj . The sample was mounted in the mixing chamber 
of a 3 He/ 4 He-dilution refrigerator with a base tempera- 
ture of 90 mK. The mobility of the cooled 2DEG was 93 
m 2 /Vs. A DC bias voltage of 10 /j,V was applied across 
the dot, and the current is measured with a resolution of 
500 fA. From capacitance measurements jj), we find the 
electronic dot area A=190 nm x 190 nm (the depletion 
length in such devices can be smaller than in structures 
defined by top gates [p0| ) . The single-electron charging 
energy is E c =1.25 meV and the spin-edgenerate level 
spacing A = 200 [ieV . 
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FIG. 1. (a) Right inset: AFM picture (taken before evap- 
oration of the top gate) of the oxide lines (bright) that define 
the dot, coupled to source (S) and drain (D) via tunnel bar- 
riers, which can be adjusted with the planar gates PCI and 
PC2. Gates I and II are used to tune the dot. Main figure: 
Conductance G as a function of Vj, showing Coulomb block- 
ade resonances. Left inset: fit (line) to one measured CB peak 
(open circles), see text, (b) Linear fit (line) of the peak spac- 
ing AV/as a function of Vi (dots). The average peak spacing 
is almost constant, indicating small shape deformations. 

The measurements have been carried out in the weak 
coupling regime, hT <§; fc^T <C A. Here, T denotes the 
coupling of the dot to source and drain. The conductance 
G was measured as a function of the voltage Vj applied 
to the planar gate I (see inset in Fig. 1(a)). Magnetic 
fields B applied perpendicular to the sample surface 
and Vjj were used as parameters. The observed CB 



oscillations (Fig. 1(a)) are fitted to a thermally broad- 
ened line shape, i.e., G(Vr) = G ma xCosh~ 2 (rjVi /2/cbT) 
yielding an electron temperature of T=120mK, as 
well as the positions and amplitudes of the peaks. Here, 
1] = O.lleV/V is the lever arm. Fig. 1 (b) shows the peak 
spacing AVj as a function of Vj. Compared to conven- 
tional dots defined by top gates, we find a much smaller 
variation of the average peak spacing as V/ is tuned, 
although the fluctuation of individual spacings is 15% of 
E c . A linear fit gives a slope of AVr/Vj = 6.7 • 10~ 4 . 
Hence, the capacitance between the dot and gate I varies 
only by 3% over the whole scan range, as compared to, 
for example, a factor of 3 in Ref. H. This indicates that 
tuning gate I or II predominantly changes the energy 
of the conduction band bottom, while the dot is only 
slightly deformed. By applying the method of Ref. |is|j 
to a hard- wall confinement, we estimate x « 0.15 for our 
dot as a lower limit. 

In Fig. 2 (a), five consecutive CB peaks are shown as a 
function of B. A pronounced pairwise correlation of both 
amplitude and peak position is observed (peak b corre- 
lates with peak c, and peak d with peak e, respectively). 
Observation of a pairing has been reported previously 
and interpreted as spin pairing Jl3] , but not been further 
investigated. 

We interpret this parametric pair correlation in terms of 
a model recently developed by Baranger et al. Q . The 
constant interaction model is used to subtract E c from 
the peak spacings. The remaining individual energy sep- 
arations equal A/2 on average and reflect the fluctuating 
level separations inside the quantum dot, which consist 
of two parts. We assume that two paired peaks belong to 
the same spatial wave function, labelled by i, of opposite 
spin, and are split by an interaction energy while 
the energy of consecutive states with different orbital 
wave functions differs by Aj — Since the separations 
between the two levels of equal spin of spin pair i and 
(i+1), Aj, and possibly also vary as a function of B, 
levels may cross and the ground state of the dot can be 
either a singlet or a triplet state. At the singlet-triplet 
transitions, kinks in the parametric peak evolution occur 
and the pair correlation is interrupted ]22| . We can iden- 
tify such kinks in our data, among other features. Fig. 2 
(b) shows the amplitudes of peaks c, d and e. The corre- 
lation between peaks d and e is very strong around B=0. 
For 0.4T<B<0.61T, this correlation is interrupted, while 
the amplitudes of peaks c and e are correlated instead. 
In this regime, correlated kinks in the evolution of peaks 
c and d are observed (Fig. 2 (c)). In Fig. 2 (d), a possible 
corresponding scenario for the parametric dependence 
of energy levels is sketched: (left) two avoided crossings 
occur between level pair i and level pair i+1 . This leads 
to the position of peaks c, d, and e as sketched in Fig. 2 
(d), right, corresponding to the difference in energy upon 
changing the electron number in the dot. Consequently, 
positions and amplitudes of peaks c and e should be cor- 
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related in 0.4T<B<0.61T, as observed. Note that this 
correlation is interrupted around B=0.5T, possibly due 
to the influence of another energy level. 
Also, S^de is not constant over the full range of B. While 
£de ~ 0.05A for B<0.22T, the positions of peaks d and 
e are not detectable in 0.22T<B<0.32T, since their am- 
plitudes vanish. As the peaks reappear, £de has jumped 
to Q e ~ 0.25A. We speculate that possibly a level cross- 
ing has occurred in the regime where the amplitudes 
are suppressed, and hence for B<0.22T, a different level 
pair is at the Fermi energy than for B>0.32T. In addi- 
tion, we note that although £ fluctuates as B is varied, a 
systematic change of £ with B is not observed, which in- 
dicates that Zeemann splitting plays a minor role. From 
the data of Fig. 2, we estimate the average interaction 
energy to £ ~ 0.5A by averaging over all peaks and mag- 
netic fields. Baranger et al. have estimated £ 0.6A 
for r s — 1. Hence, our findings can be considered as be- 
ing in agreement with existing theory, while we are not 
aware of a theoretical prediction for er^ . From the above 
phenomenology, we conclude that for dots with stronger 
shape deformations, and hence more level crossings, or in 
dots with larger r s (and thus larger £), the spin pairing 
is frequently interrupted and difficult to detect. Also, 
the Kondo effect can occur in neighboring peaks ]24| ] 
when spin pairing is interrupted. 

We proceed by discussing the effect of spin pairing on the 
NNS distributions. In Fig. 3, the measured histograms 
of the normalized NNS distributions for GOE (a) and 
GUE (b) are shown. The ensemble statistics have been 
obtained by measuring G(Vj), and either by changing 
the magnetic flux by one flux quantum cf> = — through 
the dot (GUE), or by stepping Vu in units of one CB 
period (GOE). Each individual V/-sweep contains 15 CB 
resonances in the low coupling regime. The total number 
of peak spacings used is 120 for GOE, and 210 for GUE, 
respectively. The individual level spacings s in units of 
A are obtained by using the fit of Fig. lb; its expec- 
tation value is 1=0.5. Both histograms are asymmetric 
and show no evident bimodal structure. By including 
the effect of spin pairing into the statistics, however, we 
can interpret them as modified bimodal distributions: 

(i) The (5-function in the non-interacting NNS distribu- 
tion P(s) with the expectation value of ~ss=0 (eq. (1)) is 
shifted to = £* and, as a reasonable assumption [B5| , 
broadened according to a Gaussian distribution with the 
standard deviation <j£* . Here, £* denotes the interaction 
energy in units of A. 

(ii) Since one level of a spin pair i is shifted upwards 
in energy by £j, the separation between the upper level 
of spin pair i and the lower level of pair (i+l) is given 
by Aj-fj. Consequently, P /3 (s) in eq. (1) is shifted to 
Sp/3 = l-£* and convoluted with the Gaussian distribution 
function of £*. 

Combining these two components, the modified NNS 
distribution reads 




FIG. 2. (a) Logarithmic grayscale plot of paramatric vari- 
ations in a magnetic field B for 5 consecutive CB peaks. A 
pair correlation in peak position amplitude is observed, which 
is interrupted in certain ranges of B, for example in the re- 
gion between the dashed lines, (b) Parametric amplitudes for 
peaks c,d, and e, offset by 0.2 e 2 /h each. The correlation be- 
tween peak d and e is lost in 0.4T<B<0.6T, and e correlates 
with c instead, (c) The corresponding position of the peak 
maxima. The traces are offset for clarity. At magnetic fields 
labelled by 1 and 3, kinks in the peak position occur, while 
the separation between peak d and e jumps across the region 
of suppressed amplitude from £de to £' de . (d) Scheme of a pos- 
sible double anticrossing between spin-paired level i and i+l 
(left, the black arrows indicate the spin), which could lead to 
the observed structure in the correlation for peaks c, d and e 
(right). 
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Here, the "x" denotes the convolution. Since A is deter- 
mined by the dot size and the material parameters, we 
can ntPf nt (e,^e) 
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FIG. 3. Measured NNS distributions (gray bars) for B=0 
(a) and B / (b). The bold solid curves are the fits to 
-P/nt(C*i °"£* )i with the fit results as indicated in the figure 
(see text). Also drawn are the two components of Pf nt , i.e. 
the Gaussian distribution of separations between spin pairs, 
and its convolution with the corresponding Wigner surmises. 
The inset compares the GUE data to P 2 (s) (eq. 1), using the 
spin-resolved level spacing A/2 as the average peak separa- 
tion. 



to the measured NNS distribution with the two fit pa- 
rameters £* and crj* (Fig. 3). We obtain £* = 0.65 and 
<j e = 0.35 for GOE, as well as J* = 0.53 and cr e = 0.34 
for GUE. Hence, we find that £ is higher for GOE than 
for GUE, which is in agreement the theoretical prediction 
p2| . The fluctuation of £ is found to be independent of 
the Gaussian ensemble, and does not vary continously 
with B 

In these fits, we have assumed that two electrons are 
always successively filled in one spatial wave function, 
i.e. we have neglected situations in which ^ > Aj. In- 
clusion of avoided crossings would require more clearly 
pronounced kinks than those in our data (sometimes the 
pair correlation is lost while a kink is not clearly visible). 
More experiments as well as theoretical work is necessary 
to investigate this dependence, also with respect to fluc- 
tuations in £ with B. 

Finally, we consider how our data can be modelled when 
complete absence of spin pairing is assumed. In this 
case, the mean level spacing would be A/2. Compar- 
ing a correspondingly normalized Wigner surmise to our 



data gives an extremely poor result (inset in Fig. 3): the 
measured NNS distribution appears too wide by a factor 
of w 2. 

In summary, we have observed spin pairing effects in a 
- compared to dots investigated in earlier experiments - 
rigid quantum dot with reduced electron-electron inter- 
actions. We have observed spin pairing which persists 
as a magnetic field is varied, but is interrupted by kinks 
as well as other structures in the parametric evolution 
of the Colomb blockade peaks. We have extracted the 
average interaction energy £ between states of identical 
spatial wave functions but opposite spin. Furthermore, 
we explain the measured distributions of nearest neigh- 
bor spacings as being composed of the two branches of 
a modified, bimodal Wigner-Dyson distribution, which 
takes £ and its fluctuation into account. 
It is a pleasure to thank H. U. Baranger, E. Mucciolo, K. 
Richter, and F. Simmel for stimulating conversations and 
discussions. Financial support from the Schwcizerischer 
Nationalfonds is gratefully acknowledged. 



[1] For a review, see L.P. Kouwenhoven, CM. Marcus, 
P.L. McEuen, S. Tarucha, R.M. Westervelt, and N.S. 
Wingreen, "Electron Transport in Quantum Dots", in 
Mesoscopic Electron Transport, Proceedings of a NATO 
Advanced Study Institute, edited by L.P. Kouwenhoven, 
G. Schon and L. L. Sohn, (Kluwer, Dodrecht, Nether- 
lands, 1997), ser. E, vol. 345, pp. 105-214. 

[2] M.L. Mehta, Random Matrices, Academic Press, Lon- 
don(1991). 

[3] C.W.J. Beenakker, Rev. Mod. Phys. 69, 731 (1997). 

[4] C. W. Beenakker, Phys. Rev. B 44, 1646, (1991). 

[5] M.A. Kastner, Rev. Mod. Phys. 64, 849, (1992). 

[6] U. Sivan, R. Berkovits, Y. Aloni, O. Prus, A. Auerbach, 

and G. Ben-Yoseph, Phys. Rev. Lett. 77, 1123, (1996). 
[7] F. Simmel, T. Heinzel, and D. Wharam, Europhys. Lett. 

38, 123 (1997). 
[8] SR. Patel, S.M. Cronenwett, DR. Stewart, A.G. 

Huibers, CM. Marcus, C.I. Duruoz, J.S. Harris Jr, K. 

Campman, and A.C Gossard, Phys. Rev. Lett. 80, 4522, 

(1998) . 

[9] F. Simmel, D. Abusch-Magder, D. A. Wharam, M. A. 
Kastner, and J. P. Kotthaus, Phys. Rev. B 59, R10441, 

(1999) . 

[10] SR. Patel, DR. Stewart, CM. Marcus, M. Gokcedak, 
Y. Alhassid, A.D. Stone, C.I. Duruoz, and J.S. Harris, 
Jr, Phys. Rev. Lett. 81, 5900, (1998). 

[11] R. Berkovits and B.L. Altshuler, Phys. Rev. B 55, 5297, 
(1997). 

[12] Y.M. Blanter, A.D. Mirlin, and B.A. Muzykantskii, Phys. 

Rev. Lett. 78, 2449, (1997). 
[13] A. Cohen, K. Richter, and R. Berkovits, Phys. Rev. B 

60, 2536, (1999). 



4 



[14] 
[15] 

M 

[17] 
[18] 
[19] 

[20] 



P. N. Walker, Y. Gefen, and G. Montambaux, Phys. Rev 
Lett. 82, 5329 (1999), and references therein. 



K.-H- Ahn, K. Richter, and I.-H. Lee, |cond-mat/9906046 

(unpublished), and references therein. 

Y. Alhassid, P. Jacquod, and A. Wobst, cond 



mat/9909066 (unpublished), and references therein. 
D.R. Stewart, D. Sprinzak, CM. Macus, C.I. Duruoz, 
and J.S. Harris Jr., Science 278, 1784, (1997). 
G. Hackenbroich, W.D. Weiss, and H.A. Weidenmiiller, 
Phys. Rev. Lett. 79, 127, (1997). 

R.O. Vallejos, C.H. Lewenkopf, and E.R. Mucciolo, 
Phys. Rev. Lett. 81, 677, (1998); R.O. Vallejos, C.H. 
Lewenkopf, and E.R. Mucciolo, Phys. Rev. B 60, 13682 
(1999). 

R. Held, T. Vancura, T. Heinzel, K. Ensslin, M. Hol- 
land, and W. Wegscheider, Appl. Phys. Lett. 75, 1134 



(1999); R. Held, S. Liischer, T. Heinzel, K. Ensslin, 
and W. Wegscheider, ibid. 75, 1134 (1999); S. Liischer, 
A. Fuhrer, R. Held, T. Heinzel, K. Ensslin, and W. 
Wegscheider, ibid. 75, 2452 (1999). 

[21] A simple image charge calculation reveals that the top 
gate reduces r s only by « 2%. Nevertheless, it improves 
screening of impurity potentials and the confining poten- 
tial becomes steeper. 

[22] H. U. Baran ger, D. Ullmo, L. I. Glazman, cond 



mat/9907151 



[23] D. Goldhaber-Gordon, H. Shtrikman, D. Mahalu, D. 

Abusch-Magder, U. Meirav, and M.A. Kastner, Nature 

391, 156 (1998). 
[24] S. M. Maurer, S. R. Patel, C. M. Marcus, C. I. Duruoz, 

J. S. Harris, Jr. Phys. Rev. Lett. 83, 1403, (1999). 
[25] H.U. Baranger, private communication (2000) 



5 



